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Based on the three parametric Lorenz system, a model was developed that permits to describe the behavior 
of the plasma-condensate system near phase equilibrium in a self-consistent way. Considering the effect of 
fluctuations of the growth surface temperature, the evolution equation and the corresponding Fokker-Planck 
equation were obtained. The phase diagram is built which determines the system parameters corresponding to 
the regime of the porous structure formation. 
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1. Introduction 



Nowadays, modern nanotechnologies are developed by using a variety of methods, one of which is 
the condensation process in the steady state close to phase equilibrium. This method makes it possible 
to obtain various structures of a condensate, fractal surfaces, porous structures, etc. |1, 2]. The principal 
feature of this condensation process is the plasma-condensate system being close to phase equilibrium. 
Consequently, the adsorbed atoms are arranged on the active centers of crystallization forming struc- 
tures with different architectures. In quasi-equilibrium conditions for continuous copper condensation 
of about 7 hours, the formation of highly porous structures, whiskers, and some intermediate structures 
(fibrous structures with alternating crystalline and porous parts) was observed. Of particular interest is 
the structure shown in figure[TJ which unlike ordinary single crystals, is realized under unstable temper- 
ature regime. 

These structures may be of great practical inte- 
rest. For example, they can be used as a molecular 
sieve. 

However, a question arises regarding the rea- 
sons of holding the plasma-condensate system near 
phase equilibrium. Taking into account a universal 
nature of the condensation process, we assumed 
earlier (2|-l5]] that it is caused by self-organization 
of a multi-phase plasma-condensate system. From 
the physical point of view, the mentioned self- 
organization is explained by an increase of the en- 
ergy of the adsorbed atoms which results in the 
temperature increase of the growth surface un- 
der the effect of the plasma within the condensa- 
tion process. On the other hand, an increase of the 
growth surface temperature is compensated by the 
desorption flow of the adsorbed atoms, which are 
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Figure 1. Porous structure of the copper, ob- 
tained by spray-deposited material at high dis- 
charge power ullal. 
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responsible for supersaturation. As a result, within the framework of synergetic ideology our consid- 
eration is based on the three-parameter Lorenz system |7l|8|]. 

The paper is organized as follows. In section[2] the self-organized system that forms the basis of our 
consideration is described. Section[3]discusses the statistical analysis of the motion equation. The station- 
ary solution of the Fokker-Planck equation is considered in section [4] General conclusions are presented 
in section[5] 



2. Basic equations 

Considering that a quasi-equilibrium condensation is provided on the growth surface due to a self- 
consistent development of the processes in the plasma volume, we will further use a three-dimensional 
(volume) concentration of the condensate N and a two-dimensional (surface) concentration n = Na. Here 
a is a scale factor, which plays the role of the lattice parameter and its value will be determined below. 

For a given value of the equilibrium concentration n e , the increasing supersaturation n - n e is pro- 
vided by the diffusion component defined by the Onsager relation |9, 10] for the adsorption flow 

/ ad = D|VAr|-^(Af ac -AT). (2.1) 
A 

It takes into account that the main decrease in the concentration value takes place near the cathode layer, 
whose thickness is determined by the screening length A. The latter and the diffusion coefficient D are 
given by the equations |9, 10] 

A 2 = £- , D = — , (2.2) 

where e, a are the dielectric permittivity and the conductivity of the plasma, respectively, T p is its tem- 
perature measured in the energy units; e, Ni are the charge and the total concentration of the ions of the 
deposited substance and the inert gas. 

In the second relation J2.lt . it is considered that at the upper boundary layer of the cathode the volume 
concentration of the deposited atoms is reduced to the accumulated value N ac , and the lower boundary 
of this layer presents the growth surface, near which the concentration of atoms is N. 

A decrease of supersaturation n - n e is ensured by the desorption flow J, which is directed up from 
the growth surface, so that / < 0, while the value of the adsorption flow / ac i > 0. In case there is no 
condensate (when all the adsorbed atoms have evaporated from the substrate), the condition / = -/ a( j 
is performed for the desorption component. Here, the accumulated flow / a d is defined by the equation 
d2.1> , where N - N e . The diffusion changes of the concentration N of the deposited atoms are presented 
by the continuity equation hla + VJ ad = 0. Here, the point over n denotes the differentiation with respect 
to time and the source effect is given by the estimate 

/ ad {DlA 2 ){n-n e ) 

IVJadl « -T- a " " • (2 - 3) 

A a 

Thus, the diffusion dissipation of concentration is expressed by the equation n = {D I A 2 ){n - n e ) I a. 

On the other hand, the velocity of desorption of atoms f v Ndv in volume v, based on the growth 
surface s, is as follows: 

J Ndv = -J {V])dv = -J Jds, (2.4) 

v v s 

where the first equation takes into account the continuity condition, while the second equation considers 
the Gauss theorem. As a result, the total change of concentration n = n{t) near the growth surface is 
described by the equation 

n e - n 

h=— /. (2.5) 

At the same time, the characteristic relaxation time of the supersaturation is determined by the equalities 

A 2 £ 

r n = — = - , (2.6) 

D 4na 
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the second equality being in agreement with the second relation i2.2) . 

Within the framework of the synergetic picture |8], the quasi-equilibrium condensation process is 
caused by the fact, that along with an increase of the supersaturation n-n e , the condensed atoms trans- 
fer the excess of their energy to the growth surface. As a result, its temperature T (measured from the 
ambient temperature) increases as well. This enhances the evaporation of the deposited atoms due to an 
increase of the absolute value of the desorption flow / < 0, which compensates the initial supersaturation. 

Thus, an appropriate representation of the sequential picture of quasi-equilibrium condensation pro- 
cess requires a self-consistent description of the time dependence of the concentration n{t) of adsorbed 
atoms, the growth surface temperature T(t) and the desorption flow /(f). According to @], the evolution 
equations of these values contain dissipative components and the terms presenting positive and negative 
feedbacks, the balance of which provides a self-organization process. Thereby, in the equation d2.5t . the 
first term on the right hand side represents the dissipation contribution, and the second term presents a 
linear relation between the rate of the concentration changes and the desorption flow. 

The evolution equation for the temperature of the growth surface is presented in a similar way 

T T f = -T-a T nJ + ((t), (2.7) 

where tj is a corresponding relaxation time, aj > is the coupling constant. In contrast to the equa- 
tion J2.5I . it is assumed that dissipation leads to the relaxation of the growth surface temperature to the 
value T — C0. The second term represents the nonlinear relationship of t with concentration and flow. 
Since the structure shown in figure [1] was obtained at an unstable temperature regime, the third term 
on the right hand side of <2.7) is a stochastic source of temperature changes representing the Ornstein- 
Uhlenbeck proces^]: 

I ( \t-t'\] 

<f(f)> = 0, {Ut)Ut')) = — expl-- -\. (2.8) 

Here, I is the intensity of temperature fluctuations, xr is the time of their correlation. 

To ensure self-organization, it is required to compensate the negative relationship in the expression 
<2.7t by a positive component in the evolution equation of the flow: 

t// =-(/«. +7) + a/fir, (2.9) 

where tj is a corresponding relaxation time, / ac is the accumulation flow, aj > is a constant of a positive 
feedback, allowing the growth of the / due to the mutual effect of the concentration of the adsorbed atoms 
and the growth surface temperature. 

Thus, equations <2.5t . <2.7> . <2.9> present a synergetic system, where the supersaturation n - n e is 
reduced to the order parameter, the temperature T of the growth surface - to the conjugate field, and 
the desorption flow / - to the control parameter @. As a result, the task is to investigate the possible 
stationary regimes in a stochastic plasma-condensate system, in particular, to consider the regime of the 
formation of porous structures. 

The most simple investigation of the system (2.5) , d2.7> . d2.9> is possible within a dimensionless form 
using the characteristic scales for the time f, the concentration n, the temperature of the growth surface 
T, the flow /, and for the intensity of the temperature fluctuations I: 

t s = T n , n s = a~ 2 , T s = e, J s = T~ n l a~ 2 , I^T^a] 2 , (2.10) 

where the above-mentioned length a - {ajaj) lli and energy e = {T n aj)~ l were used. 

Thus, the dimensionless system of equations describing the fluctuational transition in a plasma-con- 
densate system takes the form 

n = -(n-n e )-J, 
ef = -T-nJ + ({t), 

a J = -ihc+D + riT, (2.11) 

1 It should be noted that the condition T = does not correspond to the absolute zero since the temperature T is measured from 
the ambient temperature. 

2 Investigation of the temperature fluctuation in the form of white noise was carried out in 0|. 
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where we introduced the relations for the relaxation times 

Tt Tj 

e=—, a=—. (2.12) 

3. Statistical analysis 

While this system has no analytical solution, we will use the approximation t„ = T/ » tj, which 
means that the temperature varies most rapidly. This situation is realized in the experiment very rarely, 
but the structure, presented in figure[l] is obtained exactly under unstable temperature regime (unstable 
cooling). 

Then, on the left hand side of the equation 12. 9> , we can assume ef = 0, and the conjugate field is 
expressed by the equation T - -nj + ({t). 

After some simple mathematical operations Q[ll|[l2|] the system 12.111 reduces to the evolution equa- 
tion having a canonical form of the nonlinear stochastic Van der Pol oscillator 1 13] 

ah + Y{n)h = f(n) + g(n)((t). (3.1) 

Here, the friction coefficient j{n), the force f{ri) and the noise amplitude g{n) are presented by the 
equations 

j{n) = l + a + n 2 , 

fin) = J ac -{n-n e ){\ + n 2 ), 

g{n) = n. (3.2) 

Then, the task is to find a distribution function of the system in the phase space formed by the con- 
centration n and the rate of its change p-ah depending on time t. 

To this end, the Euler equation 13. 1\ is conveniently represented by the Hamilton formalism 

h = a~ l p, 

p = -a- l j(n)p + f(n) + g{n)at). (3.3) 

Thus, the above-mentioned probability density P[n, p, t) is reduced to the distribution function p{n, p, t) 
for the solutions of the system 13. 3> : 

P{n,p,t) = (p{n,p,t)) ( , (3.4) 

where (. . .)( means the averaging over noise (. 
We will proceed from the continuity equation 



— p{n, p,t) + \ — [np{n, p,t)]+— [pp{n, p, t)] } = 0. (3.5) 
of [on dp ' 



Further, the substitution of the equalities 13. 3t leads to the Liouville equation 



d 

p(n, p, t) = -g[n)({t)—p{n, p, t), (3.6) 
op 



where the operator 

p d d 
a an op 



(3.7) 



Turning to the interaction representation (l^l 

Q{n,p,t) = e^ {n ' p)t p{n,p,t), (3.8) 

the equation (376) takes the form 



at dp 
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where e is a dimensionless small parameter yj]. Then, using the cumulant expansion method 1 15], one 
can obtain the kinetic equation^ 



dg{n,p, t) _ 2 
dt 



t 

J {3t(n, p, t)0t{n, p, t')) {g{n, p, t'))dt', 



(3.10) 



neglecting the terms of c 3 order jg]. 

Since a physical time t is usually much longer than the noise correlation time T( , the upper limit of the 
integration can be set equal to infinity Then, returning from the interaction presentation to the original 
presentation, for the distribution function <3.4> we obtain 



— +££{n,p) 
at 



P(n, p, t) = e~ 2 J{P{n, p, t). 



(3.11) 



Here, jV is a scattering operator, which is given by the expression 



JY = [M (f)-7(ra)MiW] g 2 (n)—+cM 1 {t)g 2 (n) 



1 dg{n) j d d 2 
I h p 



g{n) dn [dp dp 2 j dndp 
where MqU) and M\{t) are the moments of the correlation function 12.81 



+ &{£ 2 ), (3.12) 



' o 



From equation 13.131 one can obtain 



(3.13) 



Molt) = I, Mi(f) = /T f . 



(3.14) 



Since, for this task, a complete distribution function P{n, p, t) has a lower practical interest than its 
integral 

&{n, f)= J P(n,p, t)dp, 
it makes sense to consider the moments of the initial distribution function 



(3.15) 



& > i[n,t) = J p'P{n, p, t)dp. 



(3.16) 



Then, the zero moment is reduced to the required integral (3.15) . 

Multiplying equation <3.1H by p l and integrating over all p, we arrive at the relation that can be 
written as a Fokker-Planck equation presented in the Kramers-Moyal form 1160 



d&>{n, t) d d 2 

— = -— [D^nWin, t)] + [D 2 {n)2?{n, t)\ , 

at an am- 



where the drift coefficient 



Diln) = 



and the diffusion coefficient 



are presented by the functions <3T2" 



g 2 {n) dy{n) dg[n) 
fW-Mo(t) ° 2 J +M 1 {t)g{n) h 



D 2 {n) = M {t) 



j 2 [n) dn 
« 2 {n) 



dn 



j 2 [n) 



(3.17) 



(3.18) 



(3.19) 



3 It takes into account that the time derivatives in equations 13.61 , 13.91 were treated according to the Stratonovich rule. 
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4. Stationary solution 



A stationary solution of the Fokker-Planck equation [lj| yields a stationary distribution (5 Hill 



&>(ri) 



D 2 {n) 



eXP i^ 



dn', 



(4.1) 



where the partition function Z is presented by the equation 



oo n 

f dn {D l {n') A , 
= / exp / -dn . 

J D 2 (n) v ] D 2 {ri) 



(4.2) 



The extremum condition for the distribution (4~1 



D l {n)-—D 2 {n) = 
on 



(4.3) 



defines the stationary states of the plasma-condensate system. 

Substituting expressions (3A8\ , (3A9i , f3TZt . and 
<3.14> into the equation <4.3> we obtain the equation 
denning the stationary concentration dependence 




Jac — 



21{l + cr)n 
[a + u) + n 2 \ 



■ + (n - w e ) (1 + n ) - lT( n. (4.4) 



Then, the condition that restricts the domain of 
the existence of the solution n — corresponding 
to the complete evaporation of the condensate from 
the growth surface, has the form 



/ac — 



(4.5) 



The corresponding phase diagram of the system 
is shown in figure [2] 

While figure [2] (a) shows the effect of the system 
parameters (the equilibrium concentration n e and 
the correlation time of fluctuations t^), figure [2] (b) 
considers in detail the domains of the phase dia- 
gram. In particular, the domain C corresponds to 
the condensation process, the domain S is charac- 
terized by the formation of porous structures, and a 
complete evaporation of the condensed matter takes 
place at the domain O. The domains, which are in- 
dicated by two letters, meet the coexistence of the 
above mentioned regimes. 

It is more convenient to understand the pro- 
cesses occurring in each domain considering the ex- 
ample of the stationary concentration dependence 
presented in figure[3] 

Each point on the phase diagram [figure [2] (b)] 
corresponds to the ray in figure[3](a), (b). 

For example, for the ray 1 [figure [3] (a)] only the 
condensation process is realized. It corresponds to 
the point C' characterized by a sufficient stationary concentration n. With a decrease of the accumulated 
flow (ray 2), there is observed a gradual disassembly of the previously formed condensate. This situation 
corresponds to the existence of two steady states with different concentrations (points C and S'). 



Figure 2. Phase diagram of the system. The solid 
line corresponds to « e = 0.25,t^ = 0.5, the dashed 
line - to n e = 0.75,Tf = 0.5, and the dotted line - 
to « e = 0.25,t> = 0.75. The letters indicate the rel- 
evant domains of the phase diagram, and the dots 
(marked by numbers) correspond to the parame- 
ters at which the stationary concentration depen- 
dence (figurefj} is analyzed. 
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Figure 3. The dependence of the stationary concentration n on the accumulated flow / ac at n e = 0.25, 
t ( = 0.5, (a) I = 8,(b)I= 14. 

It should be noted that the additional intersection point of the main curve and the ray 2 (which is 
located between C and S') applies to the non-physical plot and, therefore, it is not considerecQ. Turning 
to the state of the surface disassembly (point S'), it is worth noting that in this case the usual evaporation 
of the upper layer of the condensate does not take place. First, the atoms which are less connected with 
the crystallization centers, are detached from the condensate surface. With a further decrease of the 
accumulated flow (ray 3), the only state of the surface disassembly remains (point S). This is the situation 
that characterizes the pattern shown in figure [1] and is of great interest to us. Going to the ray 4, the 
disassembly is replaced by the usual evaporation process (point O). 

Analyzing the relationship (4.4) for higher intensity of fluctuations [figure [3] (b)], one can see that 
some changes occur. As previously, only condensation process (point C") is realized for ray 5, while ray 
6 is characterized by the coexistence of disassembly (point S) and condensation (point C') processes. The 
main difference is found for the ray 7, when together with the condensation process (point C), evapora- 
tion (point O') takes place. At the parameters specified for the ray 8, only evaporation occurs. 

5. Conclusion 

Based on the above analysis, we can conclude that processes occurring in the plasma-condensate sys- 
tem can be represented within the system <2.5I . (277) and J2.9> describing the self-consistent behavior of 
concentration, temperature of the growth surface and desorption flow. Taking into account the fluctua- 
tions of the growth surface temperature with the correlation function 12.81 makes it possible to describe 
the most specific state (disassembly of the surface), when porous nanostructures may be formed. In ad- 
dition, as shown in figure [2] (a), the system parameters have a significant effect on the domain of the 
formation of such structures. With an increase of the correlation time of fluctuations, this domain signifi- 
cantly decreases and shifts towards the lower values of the fluctuation intensity, while an increase in the 
equilibrium concentration results in a less significant decrease, as well as causes a shift along two axis 
(fluctuation intensity and accumulated flow). As a real experiment our theoretical approach has 

shown that the state of the surface disassembly is rarely realized. However, controlling the parameters 
of a system, we can reach the regime under which porous nanostructures are formed. 
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flocnifl>i<eHHfl HaHonopi/icmx MaTepia/iiB 3a yMOB 
KBa3npiBH0Barn 

O.B. rOmeHKO, T.I. >Ki/i/ieHKO 

CyMCbKuPi flepxaBHwfi ymBepcwTeT, By/i. PuMCbKoro-KopcaKOBa, 2, 40007 CyMn, YKpaiHa 

Ha ocHOBi TpnnapaMeTpw4H0i cwcreMU /lopei-ma 6y/ia po3BHHeHa MO^eyib, luo p,03B0nne caMoy3roflxeHHM 4H- 
hom oni/ica™ noBefliHKy cucreMW n/ia3Ma-KOHfleHcaT no6/in3y cf>a30B0T piBHOBam. BpaxoByK)4n Bn/inB cf>/iyKTy- 
au,\v\ TeMnepaTypn pocTOBOi noBepxHi, 6y/in 3HatifleHi piBHUHHa eBo/iKmii Ta BiflnoBiflHe piBHflHHfl <t>OKKepa- 
riyiaHKa. no6y,qoBaHa c(>a30Ba fliarpaMa, Ha ocHOBi hkoi BH3Ha4eHi napaMeipw cucTeiviw, 5=1 Ki BiflnoBiflaioTb pe- 
XMMOBi yTBopeHHfl nopucTwx CTpyKTyp. 

KnioMOBi c/ioBa: caMoopraHi3au,iR, nepecvmenm, <pa30Ba piBHOBara, KOHp,encau,\n, cpa30Ba fliarpaMa 
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